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Abstract 



The observation of a new b baryon via its strong decay into 7r+ (plus charge conju- 
gates) is reported. The measurement uses a data sample of pp collisions at = 7 TeV 
collected by the CMS experiment at the LHC, corresponding to an integrated lu- 
minosity of 5.3 fb^^. The known baryon is reconstructed via the decay chain 
— > J/i/?E^ — ?> ji^ji^ A^n^ , with A*^ — > ^pn^ . A peak is observed in the distri- 
bution of the difference between the mass of the n'^ system and the sum of the 
masses of the and 7r+, with a significance exceeding 5 standard deviations. The 
mass difference of the peak is 14.84 ± 0.74 (stat.) ± 0.28 (syst.) MeV. The new state 
most likely corresponds to the = 3/2+ companion of the S^. 
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According to the well-established quark model and corresponding spectroscopy of baryons, 
there are several predicted baryons containing one strange and one beauty valence quark. 
These include the Eb (ground state) and E^, both with total angular momentum and parity 
= 1/2+, a = 3/2+ state with angular momentum L = (often referred to, as will be done 
in this Letter, as H^), and two states with = 1/2^ and 3/2^, both with angular momentum 
L = 1. These baryons can be neutral (valence quark content u-s-b) or negatively charged (d- 
s-b). At the Tevatron, baryons with masses and decay modes consistent with the theoretical 
predictions for the ground state baryons have been observed [1 -3J, although their quantum 
numbers have not yet been established. The allowed decays of the experimentally missing Hb 
states should be analogous to the charmed sector In addition, theoretical calculations [7^ 
[Tl] predict the mass difference between the and Eb to be smaller than the mass of the pion, 
in which case the strong decay e!^ —> Eb/r is kinematically forbidden. The mass difference 
between the E^ and Eb, however, is expected to be large enough to allow such a decay. 

This Letter presents a search for the decay E^° —> Ej^7r+, with —> ]/xpE^, ]/ip —> , 
E^ — > A^n^ , and A'^ — ?> p/r^. Charge conjugate states are implied throughout. The recon- 
struction of such decays involves the presence of three secondary vertices, where the Ej^, E^, 
and A*^ decay, which are well separated from the primary interaction vertex. The analysis 
is based on a data sample of pp collisions at ^/s = 7TeV, collected in 2011 by the Compact 
Muon Solenoid (CMS) experiment at the Large Hadron Collider (LHC), corresponding to an 
integrated luminosity of 5.3 fb^^. 

The CMS apparatus is described in detail in Ref. [125. Its central feature is a superconducting 
solenoid, of 6 m internal diameter, providing a field of 3.8 T. The main subdetectors used in this 
analysis are the silicon tracker and the muon systems. The silicon tracker, composed of pixel 
and strip detector modules, is immersed in the magnetic field, and enables the measurement of 
charged particle momenta over the pseudorapidity range \rj\ < 2.5, where ?/ = — ln(tan0/2) 
and 9 is the polar angle of the track relative to the counterclockwise beam direction. Muons are 
identified in the range < 2.4 using gas-ionization detectors embedded in the steel return 
yoke of the magnet. 

The events used in this analysis were collected using the two-level trigger system of CMS. The 
first level consists of custom hardware processors and uses information from the muon systems 
to select events with two muons. The "high-level trigger " processor farm further decreases the 
event rate before data storage, requiring two opposite-sign muons compatible with being the 
decay products of a ]/tp, either promptly produced or displaced from the primary vertex (PV). 
The nonprompt J/i/? trigger requires two opposite-sign muons with an invariant mass within 
200 MeV of the ]/ip mass [13J, single muon transverse momentum pj > 3, 3.5, 4, or 5GeV, and 
dimuon vertex fit probability larger than 10% or 15%. The requirements were made tighter 
depending on the LHC instantaneous luminosity so as to limit the trigger rate. In addition, the 
dimuon vertex must be separated from the beam line by more than 3 times the uncertainty on 
the separation cTVertex/ which includes the transverse vertex resolution (independently estimated 
for each event) and the beam line position uncertainty. Finally, the requirement cos a > 0.9 is 
also applied, where a is the angle, in the plane transverse to the beam, between the dimuon 
momentum and the direction from the beam line to the dimuon vertex. The prompt J/xp trigger 
requires two opposite-sign muons with an invariant mass within 250 MeV of the ]/xp mass [13|, 
dimuon pj > 10 or 13GeV (depending on the instantaneous luminosity), dimuon rapidity 
lyl < 1.25, and dimuon vertex fit probability greater than 0.5%. 

The reconstruction of the E^^ J/'/'S^ candidates begins by identifying ]/tp — > f ^/^^ decays. 
The candidate ]/tp mesons are built by combining pairs of oppositely-charged muons (tracks 
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in the silicon tracker matched to tracks in the muon detectors) satisfying the trigger conditions 
and having an invariant mass within 150 MeV of the ]/ip mass. Candidate A*^ baryons are re- 
constructed in decays to a pion (tta) and a proton (p) with opposite charges, where the higher 
momentum track is assumed to be the proton. Both tracks are required to have a track fit 
;:\;^/ndf < 5, where ndf is the number of degrees of freedom, and at least six hits in the silicon 
tracker. The A*^ candidate vertex must have a fit x^/ ndf < 7 and be displaced from the beam 
line by more than lOcJvertex- Possible contamination from misreconstructed Kg mesons is re- 
moved by requiring that the candidate mass, when assuming that both tracks are pions, differs 
from the Kg mass |ll3| by more than 20 MeV. Candidate baryons are reconstructed by com- 
bining a A*^ candidate with a track (tts) of the same charge as the tt^. A kinematic vertex fit flU 
is performed, assuming that the track is a pion, and with the A^ candidate mass constrained 
to its world-average value ||T3|. A veto window of 20 MeV around the mass [13J is applied 
when assuming that the tts candidate is a kaon. The is then combined with the ]/ip to form 
a candidate, with a kinematic vertex fit constraining their masses to the world-average val- 
ues fl3|. On average, the events used in this analysis have around eight reconstructed primary 
vertices (PV). The PV closest in three dimensions to the trajectory is taken as the one where 
the Ej^ was produced. 

The signal selection criteria are chosen by an iterative algorithm that maximizes both the 
signal yield and the significance, calculated asS/\/S + B, where the signal (S) and background 
(B) yields are evaluated in a window of ±40 MeV around the mass, 5790.5 ± 2.7MeV [13|. 
The value of B is linearly interpolated to this signal region from the two mass sidebands, lo- 
cated between 40 and 100 MeV away from the peak value. Thirty variables are used to optimize 
the Hj^ signal selection. At every iteration, two variables are chosen randomly, where one of 
them is tightened and the other loosened. If S and the significance increase, the new set of 
values is accepted. To avoid ending up in a local maximum, an iteration is accepted if the sig- 
nificance decreases by less than a random number uniformly generated between and 10%. 
The variables used for the selection include the pj of both muons, and of the J/t/?, p, tta, tts, 
H^, and E^. The muons and can have different pj selections in the pseudorapidity regions 
< 1.2 and \ f]\ > 1.2, where \r]\ = 1.2 is the transition between the barrel and endcap detec- 
tors. The algorithm also tunes the requirements on the differences between the reconstructed 
masses and the world-average masses for the A^, S^, and ]/xp, with the latter case being differ- 
ent for \i^{]/tp)\ < 1.2 and \rj{]/tp)\ > 1.2. Also, the total J/t/; pseudorapidity coverage can be 
adjusted to be smaller than the full \ f](i/'ip) \ < 2.4 range. 

To select a sample with a good signal-over-background ratio, it is important to reject promptly 
produced tracks. This is done via a selection on their transverse impact parameter significances 
Dip/ (TDip, where the impact parameter Dip is calculated with respect to the beam line and auip 
is its uncertainty. This procedure is applied to the p, tta, and tts tracks, as well as to the A*^ and 
trajectories. Long-lived particles (A*^, H^, and S^^) are selected by requirements on their 
transverse decay lengths and on the significance Lxy /cixy of the transverse distance between 
their production and decay vertices L^y, where aixy is the uncertainty of Lxy. Minimal fit con- 
fidence levels are required on the A'^, S^, and decay vertices, while the three-dimensional 
distance significance must be smaller than the optimized thresholds for the distances between 
the trajectory and the ]/ip decay vertex , as well as between the trajectory and the chosen 
PV. 

The J/t/^H^ invariant-mass distribution for the candidates passing the selection criteria is 
shown in Fig.[l| together with the result of a fit with a Gaussian function representing the signal 
plus a second-order polynomial representing the background, which gives a signal yield of 
108 ± 14 events. The same figure also displays the invariant-mass distribution for background 
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events, in which the tth and proton have the same charge, giving further evidence that the 
observed peak corresponds to a real signal. The mass extracted from the fit is 5795.0 ± 
3.1 (stat.) MeV, in good agreement with the world-average value I13J . The corresponding mass 
resolution is 23.7 ± 3.2 (stat.) MeV, in agreement with the value 22.5 it 4.7 MeV, obtained from a 
detailed Monte Carlo (MC) simulation of the CMS detector response, using pythia 6.409 llTSi , 
EvtGen flFl, and Geant4 [IZI. 




M(J/\|/S") [GeV] 

Figure 1: Invariant-mass distributions of the selected candidates (closed circles) and of the 
J/t/^S^ pairs where the proton and tts have the same charge (open squares). The blue solid 
curve represents the result of the signal-plus-background fit to the distribution, while the 
red dashed line shows the corresponding background component. 



To search for E^" baryons, the candidates with a mass within 2.5 standard deviations of the 
fitted peak value are combined with tracks, assumed to be pions, with a charge opposite to the 
TTs charge (opposite-sign pairs) and coming from the selected PV, with a significance less than 
3 standard deviations on the distance between the track trajectory and the PV. Other quality 
requirements applied to the tracks are pj > 0.25 GeV, at least two hits in the silicon pixel layers, 
at least five hits in the entire tracker, and a track fit ;\;^/ndf < 2.5. 

The search uses the mass difference Q between the measured ]/tp 7r+ invariant mass and 
the sum of the masses of the decay products, Q = MQ/tpE^ n^) — M{]/xpE^) — M{7z), where 
M{n) is the charged-pion mass [13J. The search for new resonances in the Q distribution re- 
quires a reliable background shape. A background model is built using candidates where the 
prompt pion and the E^^ have the same charge (same-sign pairs), given that the background is 
expected to be dominated by combinatorial sources, as checked by MC studies. The measured 
momentum distributions of E^ candidates and same-sign pions (p(Eb), p{7i)), together with 
the distribution of the angle between them (a), are used to randomly generate uncorrelated 
values for p(Eb), p{n), and a. Given the limited statistical precision of the momentum 
distribution, the corresponding random numbers are generated from a parametrized version 
using the fit function f-^- (p) = p^^ e^^^P, where A:, are free parameters. The three random values 
are then combined to calculate a Q value for predicting the combinatorial background distri- 
bution. One hundred million Q values are generated in this way and the resulting distribution 
is fitted to the function Q'^^ (e^'^^Q -)_ q-<:3Q _|_ g-^Qj^ where c, are free parameters. Figure |2] (a) 
compares the Q distribution of the same-sign E^'^ candidates with the predicted background 
shape. Alternative functional forms of /n- (p) are used to estimate the systematic uncertainty 
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associated with this method, which contributes to the determination of the background param- 
eters. 

The measured opposite-sign Q distribution is displayed in Fig. |2] (b) for the range 0-50 MeV. 
The 21 events observed in the region 12 < Q < 18 MeV represent a clear excess with re- 
spect to the expected background yield of 3.0 ± 1.4 events, evaluated by integrating the back- 
ground function in this Q window. An unbinned maximum-likelihood fit is performed to 
the opposite-sign Q distribution with a Breit-Wigner distribution convolved with a Gaussian 
function, added to the background function previously described. The Gaussian resolution 
of the peak is constrained to 1.91 ± 0.11 MeV, as determined in the signal MC simulation, 
and the background parameters are allowed to float within their total uncertainties (statisti- 
cal plus systematic, added in quadrature). Figure |2] (b) also shows the result of the fit. A 
peak is clearly visible above the background continuum. The fitted parameters of the peak 
are Q = 14.84 ± 0.74 (stat.) MeV and Breit-Wigner width F = 2.1 ± 1.7 (stat.) MeV. The fitted 
Breit-Wigner width agrees with F = 0.51 ± 0.16 MeV, the value obtained following Eq. (102) of 
Ref. |18| , based on lattice quantum chromodynamic calculations. 

To evaluate the significance of the signal, the likelihood Ls_|_(, of the signal-plus-background 
fit is determined. The fit is then repeated using the background-only model to obtain a new 
likelihood Lj,. The parameters obtained from the background fit are allowed to float so that 
their uncertainties and correlations contribute to the calculation of the significance. The loga- 
rithmic likelihood ratio A/ln(Ls_|_fc/Lj,) would indicate a statistical significance of 6.9 standard 
deviations (a), corresponding to a probability of 2.5 x 10^^^ for a background fluctuation of 
this significance or more to be observed. The significance remains the same if the fit is repeated 
using the theoretically expected width, allowing it to vary within the range of the theoreti- 
cal uncertainty. The signal significance is also evaluated by generating pseudoexperiments in 
which the background distribution is varied within its statistical uncertainty, and determining 
the background fluctuation probability ("p-value") as the number of experiments that give a 
fit with the same significance or higher than in the data. The "look-elsewhere effect" |19| is as- 
sessed by searching for a peak, of width F between and 25 MeV, in the extended mass range 
< Q < 50 MeV, where the E^^ is theoretically expected. The resulting background fluctuation 
probability is p = 1.3 x 10^^, which corresponds to a 5.7(7 equivalent Gaussian significance. If 
the search range is further extended to < Q < 400 MeV, the equivalent Gaussian significance 
becomes 5.3(7". 

This analysis has been repeated using simulated B+, B'^, Bg, and Ab samples, obtained by the 
detailed MC simulation of the CMS detector already mentioned. The samples contain events 
in which the b hadron is forced to decay to a J/t/j, which decays to No evidence of peaks 

due to partially reconstructed b-hadron decays is observed in these samples. The opposite-sign 
Q distribution obtained using candidates from the lower- and higher-mass sidebands of the 
signal peak also shows no excess, indicating that the observed peak is not caused by fake 
candidates. 

The systematic uncertainty on the measured Q value is evaluated through a detailed MC sim- 
ulation. The reconstructed Q value in MC simulations is measured to be 0.23 ± 0.10 MeV above 
the generated value. This is consistent with the observation that the measured A*^ and 
masses are 0.16 ± 0.05 and 0.18 ± 0.14 MeV, respectively, above their world averages. The sum 
in quadrature of the shift and its statistical uncertainty, 0.25 MeV, is considered as the system- 
atic uncertainty due to this effect. As an extreme fitting scenario, a flat function is used for 
the background shape, leading to a Q value 0.12 MeV higher than the value measured with 
the nominal background model. Adding in quadrature this uncertainty with the previous one 
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Figure 2: (a) Same-sign Q distribution (closed circles) and result of a fit with the background 
model (red dashed curve) in the range < Q < 400 MeV. (b) Opposite-sign Q distribu- 
tion (closed circles) in the < Q < 50 MeV range, along with the result of the signal-plus- 
background fit (blue solid curve); the background term is also shown (red dashed curve). 



results in a total Q systematic uncertainty of 0.28 MeV. 

The observation of this resonance, corresponding to the one observed in the charm sector lU, 
and its mass measurement add valuable information to the understanding of the interactions 
between quarks within a baryon. 

In summary, a new Eb baryon has been observed in pp collisions at = 7 TeV, using data 
collected by the CMS experiment, corresponding to an integrated luminosity of 5.3 fb^^. The 
signal is observed with a significance exceeding 5 standard deviations. The measured Q = 
M{]/xpZ-Ti+) - M{]/xpZ-) - M{n) value is 14.84 ± 0.74(stat.) ± 0.28 (syst.) MeV. Given the 
charged-pion and masses [13], the resulting b-baryon mass is 5945.0 ± 0.7 (stat.) ± 0.3 (syst.) ± 
2.7 (PDG) MeV, where the last uncertainty reflects the present accuracy of the E^^ mass from the 
Particle Data Group |13|. While the width of the new baryon is not measured with good statis- 
tical precision, it is compatible with theoretical expectations IITSl . Given its measured mass and 
decay mode, the new baryon is likely to be the H^^, with = 3/2+. 
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